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We calculate the energy eigenvalues and eigenstates corresponding to coherent single and multiple 
excitations of an array of N identical qubits or two-level atoms (TLA's) arranged on the vertices 
of a regular polygon. We assume only that the coupling occurs via an exchange interaction which 
, depends on the separation between the qubits. We include the interactions between all pairs of 

qubits, and our results are valid for arbitrary distances relative to the radiation wavelength. To 
illustrate the usefulness of these states , we plot the distance dependence of the decay rates of the 
£SJ , n = 2 eigenstates of an array of 4 qubits, and tabulate the biexciton eigenvalues and eigenstates, and 

absorption frequencies, line widths, and relative intensities for polygons consisting of N = 2, ■ ■ • ,9 
TLA's, in the long-wavelength limit. 
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In this paper, we calculate the eigenstates of N identical (coherently excited) two-level quantum systems arranged 
■ on the vertices of a regular polygon. Such systems are known as qubits to the quantum information community and as 
two-level atoms (TLAs) to the quantum optics/spectroscopy community. The coherent excitation of identical TLAs 
[~^. , has long been of interest to spectroscopists, in connection with the theory of molecular excitons |jj for example, or 
• the phenomenon of superradiance [El. More recently, the interest has been in connection with the optical properties 
of molecular clusters or aggregates (p) , many of which properties are believed to be related to the coherent interaction 
of the aggregates with the radiation field. At the same time, multiparticle entangled states of qubits have become 
an active area of study in the field of quantum information theory. The results presented here are of relevance to 
this community in the studies of decoherence-free subspaces [|] and investigations into the entanglement properties 
of rings of qubits || . 

We emphasize the complete generality of the majority of the results obtained herein: our results are applicable to 
all systems in which excitation is exchanged between the pairs of interacting qubits. Such exchange interactions occur 
widely: For example, our theory is applicable to systems in which the coupling is via a spin-exchange interaction, or 
via a retarded dipole-dipole (quadrupole-quadrupole) interaction, such as exists in coherent dipole Q (quadrupole 
0) radiative excitation of atoms or molecules. We do not make the common approximation of including only nearest- 
neighbour interactions, but rather we diagonalize the full Hamiltonian, and for arbitrary distances relative to the 
, radiation wavelength. This is important for many physically realistic systems, in which the coupling between non- 
k> nearest neighbours can exceed that between adjacent qubits. 

The eigenstate calculations are presented in section |. In section [ A , we begin by reviewing the calculation of the 
eigenstates for single (n = 1) excitations of a system of N qubits arranged at the vertices of a regular polygon and 



interacting via an exchange interaction, valid for arbitrary N. Next, in section [B we present a method for calculating 
the eigenstates for double (n — 2) excitations of the system, also for arbitrary N. Finally in section IC, we outline 
the calculation of the triplet (n = 3) eigenstates for N — 6 and 7, and present in Tables I-III the complete set of 
eigenstates for all regular polygons up to and including N = 6; results for N — 7 are available upon request. 

In section |0|, we present some results specific to the physical realization of qubits in terms of TLAs interacting via 
a retarded dipole-dipole (or quadrupole-quadrupole) interaction. We have special interest in the total decay rates 
of these eigenstates in order to identify particularly long lived states, which may be useful in encoding quantum 
information. To quantum information theorists, these are known as "decoherence-free" states, and to spectroscopists 
as "subradiant" states ||. In general, complete subradiance exists only in the small sample limit, when distance 
effects are ignored. Since our calculations contain the complete distance dependence, they can be used to examine 
deviations from the "long wavelength" or "equal collective decoherence" assumption commonly made in the theory 
of decoherence-free subspaces. 
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In the spectroscopy community, the study of collective atomic phenomena is many years old, beginning with Dicke's 
pioneering article |gj ; for the early work, see [§, , and references therein. A detailed study of the cooperative emission 
by a fully-excited system of 3 identical atoms in some specific geometrical configurations was performed by Richter 
while the complete eigenstates for two- and three-atom systems of arbitrary geometrical arrangement can be 



found in reference 10 . The single-excitation eigenstates of linear chains were presented in of two-dimensional 
arrays in fl^| , and of rings and regular polygons in Single and double excitations of regular polygons in the 

long- wavelength limit were considered by Spano and Mukamel ||; however, they included in their calculations only 
nearest-neighbour interactions, so that our energy eigenvalues and eigenstates differ considerably from theirs. 



I. THE CALCULATIONS 



We consider systems of N identical qubits located at positions rj, each with ground state \0i), excited state 
and transition frequency uj. The free Hamiltonian is given by 

iV 
i=l 

Henceforth, we will label states in the "computational" basis, i.e. the bare uninteracting states, according to which 
atoms are excited therein. For example, the state of the N = 5 system in which atoms 2 and 5 are excited is written 
by quantum opticians as \gegge), by quantum information theorists as |01001), and by us here as |25). The state with 
all qubits in the state |0) (|1)) is denoted by \G) (\E)). 

The generic (excitation-) exchange interaction Hamiltonian of the qubits is given by 

JV 

Htm = ]T i>n,,s;s i . (i) 

i,J = l 

where Sf and S~ are the raising and lowering operators of qubit i. The sole assumption we make regarding the 
interaction potential fiy is that it is a function only of the separation between qubits i and j, = — r,-. We 
focus in this paper on qubits arranged at the vertices of regular polygons, and number them sequentially around the 
polygon (see Fig. 1). For nearest-neighbour qubits, we define f2j,i±i = a; similarly, Slis±2 = b; for N atoms there are 
L-/V/2J characteristic interactions, which we label sequentially alphabetically. 
In analysing the system of interacting qubits we are faced with two possibilities; in this paper we follow option (ii): 

(i) We can take the Oy to be real, for example equal to the well known expression for the dipole-dipole interaction 
energy. Diagonalizing the interaction Hamiltonian yields eigenvalues which are the energy level shifts. The dynamics 
of the system can then be analyzed using a master equation, which would include terms containing the dipole-dipole 
damping. By solving the master equation we can calculate various quantities of interest, in particular the total decay 
rate from a given eigenstate to all the states below it; this is useful for identifying long lived states desirable for 
quantum computing. 

(ii) Alternatively, we can include the free-atom radiation damping Hd = X^iLi ^-T 1 1 i ) { 1 i U an d as well take the fly to 
be complex. The real part of Clij is then the interaction energy, the imaginary part the inter-qubit damping. Although 
this results in a non-Hermitian Hamiltonian, it has the advantage that the imaginary part of each resulting eigenvalue 
automatically contains the total decay rate for that eigenstate (the real part is still the energy shift). This is proven 
elsewhere @, |§. 

The full Hamiltonian to be diagonalized is represented by a 2 N x 2 N matrix. Fortunately, it is block-diagonal in 
structure, breaking up into a series of submatrices, in each of which the coupled subsets of states all have the same 
number n of excited qubits. The submatrices are of dimension (^) , and the submatrix for n excited qubits is the 
same as that for N — n excited, a general property of exchange interactions; this halves the amount of work we must 
do (and we consequently tabulate results only for n = 1,. .. , L^J)- The n = and n = N eigenstates are just |G) 
and \E) respectively. 



A. n — 1 : Single excitation eigenstates 

The single excitation (or n = 1) eigenstates of a system of N qubits arranged at the vertices of a regular polygon 
were calculated years ago |l3] , guided by the symmetry of the system under rotation about an axis perpendicular to 



3 



the polygon plane; the n = 1 eigenvalues and eigenstates for N = 1 — 6 are there tabulated. Here we rewrite these 
calculations in a notation which allows us to extend them to states containing higher numbers of excited qubits, using 
the case of N = 5 as an example. 

In the subspace spanned by the basis vectors {|1), |2), • • • , |5)}, the matrix to be diagonalizcd has the form 



H (1) = 



( a b b 

a a b 

b a a 

b b 

\a b b a / 



b 

a a 



We introduce the matrix P, a generator of the 5-dimensional representation of C5 (the cyclic group of order 5): 



P = 



/o 1 \ 

10 

1 

1 

V 1 / 



(2) 



The eigenvalue equation of P is given by 



Pu 



(V) 



X v u 



(«). 



where A = e 5 



= {X v , X 2v , X 3v , X Av , X 5v ), and v = 1, • • • , 5. We define the polynomial M{x) = a(x + x 4 ) + 



b(x 2 + x 3 ), in terms of which H^l = M(P). Since H-^' t is a sum of powers of P, the eigenvectors U( v ) of P will be 
eigenvectors of as well, and we write the eigenvalue equation 



H int U (v) 



: m,„>u 



where the eigenvalues — M(X V ) = G^j + iF^. There is 1 non-degenerate eigenvalue corresponding to v = 5, 
771(5) = 2a + 26, and (5-l)/2 degenerate pairs of eigenvalues, corresponding to roots which are complex conjugates of 
each other: X v = (A 5_t) )*. The eigenvector corresponding to 777,(5) is simply ut^ = (1, 1, 1, 1, 1). For the eigenvectors 
corresponding to the degenerate pairs of eigenvalues, we choose the real linear combinations of ui v ) and W(5_„), 



RU, 



(v) 



IU, 



(v) 



-{u [v) +U( 5 -„)) 
— (tt(„) - U( 5 _„)). 



(3) 
(4) 



Together with U( 5 ), these form an orthogonal basis set for the 77 = 1 subspace. They are listed in Table II. 



B. n — 2: Double excitation eigenstates 

1. Odd values of N 

We continue with the example of N — 5 to demonstrate how to calculate the 77 = 2 (biexciton) eigenstates for 
general odd values of N. The subspace corresponding to N = 5, 77 = 2 has 10 basis states, which we take in the order 
{|12},|23},--- , |51>; |13>, 1 24> , - - - ,|52)}. 

If we define the four polynomials M u (x) = b(x+x i ) 1 Mi 2 (x) = a(x i +x 5 )+b(x+x 3 ) 1 M 2 i(x) = a(x+x 5 ) + b(x 2 +x 4 ) 
and M 2 2(x) — a(x 2 + x 3 ), then the interaction can be represented by the 10 x 10 matrix, 

- M(P\= ( M ^ P ) M 12(P)\ 

(2) 

Thus, Hl nt is partitioned into a 2 x 2 array of square submatriccs, each of dimension 5x5. The ability to write the 

matrix in this form is directly due to the ordering of the basis vectors, which allows the rotational symmetry of the 

(2) 

pentagon to be reflected in each of the submatrices. It is easy to show that for any odd value of N, H- n ' t can be 
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partitioned in this way into an array of (N — l)/2 x (A — l)/2 square submatrices, each of dimension A x A. This 
results in a dramatic simplification of the problem, for instance here we need diagonalize only a 2-dimensional matrix 
instead of the original 10-dimensional one. 

As with the n = 1 case discussed above, each matrix Mij(P) is a linear combination of P and its powers, and 
therefore has the eigenvalue equation 

Mij(P)u( v) = Mij{\ v )u [v) , 

where X v and ur v \ are the eigenvalues and eigenvectors of P. In order to obtain the eigenvalues and eigenvectors of 

(2) 

HlJ t , we first solve the eigenvalue equation 

M{x)V{x) = n(x)V(x), 

where V(x) is an eigenvector and fi(x) an eigenvalue of the two-dimensional matrix M(x). The solutions are easily 
found to be 

^(s) = 5 l M n(x) + M 22 (x) ± R[x)] , 

where 

R{x) = yj (Mnfx) - M 22 (x)f + AM l2 {x)M 21 (x), 

and 

± _ / M u (x) - M 22 (x) ± R{x) \ 
V [x >-\ 2M 21 (x) ) ■ 

The eigenvalues and eigenvectors of HjQ can then be shown by direct substitution to be {/i ± (A")} and 

U {v) — V (A )®u {v) - {y±( X v )u{v) J , 

where v = 1, • • • , 5. As with n — 1, for degenerate eigenvalues we form the real linear combinations of the eigenvectors; 
the complete orthogonal basis set is listed in Table II. 

In general, the eigenvalue equation for the n — 2 excitations of any odd-Af array of qubits is solved in the same 
way: 

(2) 

(i) The interaction matrix H\ r ' t is partitioned into an array of square submatrices, each of dimension N x N. 

(ii) The eigenvalue equation of matrix M^- (P) is solved, where P is the N x N matrix analogous to Eq. (||) . 

(iii) The eigenvalue equation is solved for the corresponding (A — l)/2 x (N — l)/2 matrix M(x), yielding eigenvalues 
{^(i)(x)} and eigenvectors {Vu)(x)}. 

(2) (2) (2) 

(iv) The eigenvalues and eigenvectors of H- n ' t are then given by (A 11 ) = G\L +iF>L} and {U^ = (\ V )®U( V )}, 

where v = 1, • • • , N, i = 1, • • • , \_(N — 1)/2J, A = e 2 ^, and the vectors {ut v ^} are the eigenvectors of the matrix P 
corresponding to the A-sided polygon. 



2. Even values of N 

The calculations for the n = 2 energies and eigenstates for even values of A cannot be described (or performed) so 

('2') 

succinctly. This is due to the fact that (A — l)/2 is an odd half-integer. As a result, the matrix for H\^ t consists of 
2 parts: an inner "core" of [(A — 1)/2J x [(A — 1 )/2j square submatrices, each of dimension A x A, plus an outer 
section of N/2 extra columns to the right and rows at the bottom of the core. For example, the n = 2 interaction 
matrix for A — 4 is given by 





b 
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a a\ 


b 





b 





a a 
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a a 
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a a 
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with an inner core matrix M(P) = b(P + P 3 ), where P is now the 4-dimensional analogue of Eq.(||). 

The calculations are performed in the following manner; we illustrate the general procedure with the example of 
N = 4: 

(i) The energy eigenvalues and vectors of the "core" matrix are obtained, in exactly the same way as described in the 
previous section for odd values of N. 

For the case of N = 4, the eigenvectors of M(P) are the same as those of N = 4, n = 1, which in turn are the same 
as those of P. They appear in Table 

(ii) These eigenvectors are then divided into 2 groups, according to their symmetry or antisymmetry. The vectors 
(1, 1, • • ■ ,1,1) and (1,-1,1,-1, ■•■,1,-1) are always eigenvectors, the former symmetric, the latter antisymmetric; 
the remainder are classified according to their symmetry under rotations of tt about the symmetry axis. 

In the case of N = 4, three of these eigenvectors (those corresponding to v = 1, 2, and 3 as listed in Table |) are 
antisymmetric, while that corresponding to v = 4 is symmetric. 

(iii) The antisymmetric eigenvectors are appended with N/2 O's; the resulting vectors are eigenvectors of H^ n L and 
the corresponding energies are found by direct substitution. 

In the case of N = 4, by appending two O's to the ends of the antisymmetric vectors, we obtain the following three 
eigenvectors of H 



(2). 



/ 1 \ 



U(i) 



1 



-1 




V o / 



(2) 



1 

-1 
1 



V o / 



( \ 

1 



-1 



V o J 



The corresponding eigenvalues are found by substitution. By symmetry, we see that a fourth (antisymmetric) eigen- 
vector of H^ t isU= (0, 0, 0, 0, 1, -1). 

(iv) The symmetric eigenvectors are extended into the rest of the n = 2 subspace in a symmetric fashion. 

For the example of N = 4, the remaining 2 eigenvectors are found from the symmetric eigenvector U( 4 ) = (1, 1, 1, 1) 

(1, 1, 1, 1, x, x), obtaining 2 (inde- 

b±R, 

They are listed together with 



of H^ t . We substitute into the eigenvalue equation for H^ t the trial vector U 

pendent) equations for x and the eigenvalues //: 26 + 2ax — /i and 4a = [ix. These have the solutions fi 
l ^£H > where R = Vb 2 + 8a 2 . This completes our set of 6 eigenvectors of H^ t 



x 



their corresponding eigenvalues in Table |I| 

We point out that the N = 4, n — 2 eigenstates are the first which depend on the actual strength of the interaction, 
and not merely on its symmetry. In Fig. 2, we illustrate the distance dependence of their decay rates. In the long- 
wavelength limit , three of the eigenstates have their decay rates unchanged from the noninteracting value of 27, and 
one state is superradiant, with an asymptotic value of 5.9307 ( see Table V). The remaining two states are subradiant: 
One shows weak optical activity, with a limiting decay rate of 0.0707, and one is completely subradiant, with decay 
rate — ► 0; thus, this state is of possible interest for the encoding of quantum information. (We have cut off the figure 
at X/r — 10 in order to retain the visibility of some of the oscillations at low values of the argument, corresponding 
to shorter wavelengths.) 



C. Triple excitation eigenstates 

To complete the sets of eigenstates for the N = 6 and N — 7 polygons, we require those corresponding to the n = 3 
excitations. These are obtained with methods very similar to those used for the n = 2 states. For N — 7, the IfS 
matrix is first partitioned into a 5 x 5 array of square submatrices, each of dimension 7x7. The solution then requires 
the (preliminary) diagonalization of a 5 x 5 matrix M(P), but is otherwise a direct extension of the method used for 
ri = 2. 

For N = 6 we choose the basis vectors in the order: {|123), |234), • • • , |612); |124), |235), • • • , |613); 

1 134), 1 245), • • • , 1 623); 1 135) , |246)}. Doing so we find that the matrix H^ t consists of a core array of 3 x 3 sub- 
matrices, each of dimension 6x6, together with an outer section of 2 columns to the right and 2 rows at the bottom 
of the core. The solution consists of 2 stages: In the first stage, the eigenvalues and eigenvectors of the core matrix 
are found, and in the second the symmetric and antisymmetric eigenvectors of the core are extended to become those 
of the full 20 x 20 matrix, in a manner entirely analogous to that employed for N = 6, n = 2. The complete set of 
eigenvalues and eigenvectors for N=6 is listed in Table III. Those for N=7 are available upon request. 
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II. COOPERATIVE RADIATIVE TRANSITIONS 



In this section, we focus on systems of identical two-level atoms (or molecular monomers) undergoing cooperative 
radiative transitions; the interatomic potential which applies in this case is the retarded multipole-multipole interaction 
EJ , The strongest and most common of these are 1-photon transitions due to the electric dipole moment operator; 
however, the same analysis can also be made for magnetic dipole or (with different f2^) higher electric multipole 
transitions |l5) as well, and even for 2-photon transitions Jl(| . We point out that in many systems an interaction 
exists between nearest neighbours (e.g. due to atomic overlap) in addition to the electromagnetic exchange interaction 
which occurs between all pairs. These forces can be included trivially in our analysis, simply by incorporating them 
into the nearest-neighbour interaction fi^jii = a. 

For the simple case of linear transition dipoles, of transition strength /i, oriented (parallel to each other and) 
perpendicular to the plane of the ring, flij can be written in the form 



Slij = i-f 



~4 2) (^)+4 2) (^«) 



where hffl is a spherical Hankel function of the second kind |17j and 7 is half the atomic Einstein A coefficient, 

2HV 5 
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For linear transition quadrupoles oriented perpendicular to the plane of the ring, the interaction is 



n 



ij — l lq 



9 



,(2), 



f 



(2), 



where "f q is half the Einstein A coefficient for the (quadrupole) transition, 



\q\ 2 m 5 
I5hc 5 



jj], However, the analysis is also valid for any system in which the transition moments of all identical units are 
oriented symmetrically, i.e. they form the same angle with the ring. For example, there exist molecular aggregates 
in biology known as "light-harvesting complexes" , in which large identical building blocks or monomers are arranged 
symmetrically in rings with an N-fold symmetry axis, whose electronic excitations have been found to extend coherently 
over the entire ring [fl8|| . The direction of the transition moment of each individual monomer forms the (same small) 
angle 8 7^ 90° with the tangent to the ring at its position. In this case, the dipole-dipole interaction is (slightly) more 
complicated in form, and is given by 



Qij = i-y 



{3(Ai • fii)(Aj ' hj) - h ■ tJ>j}hf\krij) + (fk ■ fij) h { 2> (kr i:j ) 



.(2), 



but all our calculations remain valid. 

The dynamics of the system are governed by the Lehmberg-Agarwal master equation Jl9| , which may be solved by 
projection onto any complete set of basis vectors. However, the "natural" set for this projection are the eigenstates 
of the Hamiltonian H = Hf + Hd + Hint- As demonstrated previously for N — 2 and 3 JhJ and elsewhere for general 
N pp| ], these states have the following properties: 

1. The real part of the (complex) eigenvalue gives the shift of energy of the exciton due to local field effects. 

2. The imaginary part of the eigenvalue gives the total decay rate or inverse lifetime of the state; this total decay 
rate is the sum of the individual decay rates to all states in the energy manifold below, as calculated using the master 
equation and in agreement with the total energy radiated by the system in a transition between the two states j2(J . 

3. The eigenstates form the basis set within which the population dynamics and spectroscopic properties of the system 
are most conveniently studied. 

A calculation of the n = 1 and n — 2 eigenstates of regular polygon systems of odd N for electric dipole interactions 
in the long- wavelength limit has in fact been carried out Q , but in that reference the authors included only nearest- 
neighbour coupling for the real part of the interaction. For the n = 1 subspace, the resulting eigenstates are the 
same as ours (which however include interactions between all neighbours and are valid for arbitrary wavelength), 
but the energies are very different: This difference is illustrated in Table IV, where we list the n = 1 energies in the 
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long- wavelength limit for polygons having N — 5 and N — 6, for nearest-neighbour interactions only, for linear dipoles 
with all neighbours included, and for linear quadrupoles with all neighbours included. (All energies are expressed 
in units of the static interaction energy Vn between a pair of nearest neighbours.) For the n = 2 subspace, the 
eigenstates themselves are very different from those obtained when only nearest-neighbour interactions are included, 
and a numerical comparison of the energies is not meaningful. Because the retarded interactions are intrinsically 
long-ranged, a correct calculation of the eigenstates of the physical system must include interactions between all atom 
pairs: Only in these states will the "local field" shifts be the same, and only in these states will the damping be 
the same for all atoms, so that no dephasing occurs during the evolution in time f2l|| . As well, in some systems the 
energy of interaction between second (or higher) nearest neighbours can actually exceed that between adjacent pairs 
(depending on the relative phases of the moments in the given eigenstate, and/or on the relative orientations of the 
transition moments and r^). 



The detailed emissive properties of these systems will be presented elsewhere 20 ; however, some simple properties 



are immediately evident in the eigenvalues and eigenvectors. For example, in the long wavelength limit only one n = 1 
state is optically active in absorption and emission, and it is superradiant, having an eigenvalue whose imaginary part 
— > A7; the N — 1 other single-excitation states are subradiant, with the imaginary parts of their eigenvalues — > 0. In 
general, the n = 2 states decay into n — 1 states (although these in turn may be subradiant); however, for even-iV 
polygons, there is (at least) one n = 2 state which is itself completely subradiant. 

As a simple illustration of the emissive pattern, in Fig. 2 we display a complete energy level diagram for N = 3 
|p2| . All atomic separations are equal, so that the system is described by a single interaction potential a = r y(g + if). 
On each state is indicated its total decay rate, while individual decay rates between states are indicated on the dashed 
lines connecting them. For example, the n — 2 eigenstate (1, 1, 1) has energy h(2u; + 2-fg) and a total decay rate of 
27(1 + /); it decays at a rate of 7(4 + 8/)/3 to the (symmetric) (1, 1, 1) state, and at a rate of 7(1 — /)/3 to each 
of the 2 antisymmetric states in the n = 1 manifold below it. Similarly, each of the other two n = 2 eigenstates has 
energy H(2lj — jg) and a total decay rate of 7(2 — /). The small sample/long wavelength limit corresponds to / — > 1. 
In this limit, all the n = 2 eigenstates decay, but the two antisymmetric n = 1 states are (completely) subradiant. 



A. Absorption from an external field 

If a system in its ground state is placed in a weak external field of wave vector k and polarization e\, only the n = 1 
states are excited, with a relative probability proportional to |(it(„)| J2i ^tP' ®a e ik Ri |G) | 2 . If the field is sufficiently 
intense and the losses sufficiently low, population can remain in the {ur v -\} states for long enough to allow excitation 
of the n = 2 states; and so on. 

In recent years, there has been interest in the excitation of the exciton and biexciton states of the light-harvesting 
complexes, in connection with the calculation of their third-order nonlinear optical susceptibilities The complexes 
discovered so far have diameters of the order of 10 nm, and their absorption frequencies correspond typically to 
wavelengths > 400 nm, so that the long-wavelength limit applies. In this limit, the dependence on k in the absorption 
probability is negligible, and it is easy to verify that only those states {«(„)} which are totally symmetric in the atomic 
positions are optically active, namely the states it(jv) = (!)• *• ; !)■ This gives rise to lorentzian (exciton) absorption 
lines, centred at the shifted frequencies G?^\ = ^0 + with (natural) widths 2i^l — 2Nj. 

(2) (2) (2) 

We denote the energy of state U( v ^ by G)JL = 2wo + AG, vi -., and its decay constant by F}L. The excitation of 
state U( V i) from a ring in the state U(jv) then occurs at frequency ~ Gvjyy It can be shown |2(| that the width 

of the absorption line is 2(iv^l + F^), and that in the long- wavelength limit the relative intensities of the biexciton 

(2) 

absorption lines are simply given 

In Table V we list the n — 1 shifts AG^, the (unnormalized) n — 2 eigenvectors {U^}, and the n = 2 shifts 

(2) (2) 

AG jj and decay constants F^l) ^ or ^ — 2, • • • , 9, in the long-wavelength limit. The vectors {U^} correspond to 
basis states arranged in the order {|12), |23), ■ • • , \N1); |13), |24), • • • , \N2); etc}. In Table VI we list the corresponding 
biexciton excitation frequency shifts, (natural) half widths, and relative intensities. All frequencies are expressed in 
units of Vjv, and widths in units of 7. 
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III. CONCLUSIONS 

We have calculated the eigenstates corresponding to coherent single and multiple excitations of an array of N 
identical TLA's or qubits arranged on the vertices of a regular polygon. Their coupling occurs via an exchange 
interaction which depends only on the separation between the qubits. We include the interactions between all pairs, 
and our results are valid for arbitrary distances relative to the radiation wavelength. To illustrate the usefulness 
of these eigenstates, we plot the distance dependence of the decay rates of the n = 2 eigenstates of a system of 4 
qubits arranged on the vertices of a square, and tabulate the biexciton eigenstates and eigenvalues, and absorption 
frequencies, line widths, and relative intensities for polygons of N = 2, • ■ • ,9 identical TLA's, in the long-wavelength 
limit. The states will be used elsewhere to study the emissive properties of these systems and to calculate the 
amount and distribution of entanglement in these "natural" entangled states at both zero and finite temperature. 
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FIG. 2: Distance dependence of the decay rates of the N = 4, n — 2 eigenstates 



N 


n 


V 


Eigenvalues 


Eigenvectors 


2 


1 


1 


—a 




(-1,1) 






2 


a 




(1,1) 


3 


1 


1,2 


— a 




(-1,-1, 2); (1,-1,0) 






3 


2a 




(1,1,1) 


4 


1 


1,3 


-6 




(0,-1,0,1); (1,0, -1,0) 






2 


-2a + 6 




(-1,1,-1,1) 






4 


2a + b 




(1,1,1,1) 




2 


1,3 





(o,- 


-1, 0, 1, 0, 0); (1, 0, -1, 0, 0, 0); (0, 0, 0, 0, 1, -1) 






2 


-26 




(-1,1,-1,1,0,0) 






4 






(l,l,l,l,a; ± ,x ± ) 










R = 


V6 2 + 8a 2 x± = (-b±R)/2a 



TABLE I: Eigenvalues and eigenvectors of a single pair (N=2), equilateral triangle (N=3), and square array (N=4) of qubits. 
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FIG. 3: Energy levels and decay rates for an equilateral triangle of qubits (in units of the single qubit decay rate); g (/) is the 
real (imaginary) part of the inter-qubit potential a. 



n 


V 


Eigenvalues 


Eigenvectors 


1 


1,4 
2,3 
5 


2(c2a — Cib) 
2(-cia + c 2 6) 
2a + 26 


(C2, C4, C6, Cs, Clo); (S2, S4, Sg, Sg, S10) 
(C4, Cs, Cl2, Ci 6 , C20); («4, S8, Sl2, S16, S20) 
(1,1,1,1,1) 


2 


1,4 
2,3 
5 


E±(-aa, C2b) 

E±(c2a, —cib) 

E+(a,b) 
E-(a,b) 


(C2, — Cl, — Cl, C2, 1; V±C2,V±,V±C2, —V±Cl,—V±Cl) 

(-1, -2c 2 , 2c 2 , 1, 0;v±, 0, -«±, -2w±c 2 , 2t;±c 2 ) 

( — Cl, C2, C2, — Cl, 1; — W±C1, W±, —W±Cl,W±C2,W±C2) 

(— 2c2, 1, —1, 2c2, 0; 2w±C2, 0, — 2iu±C2, w±, — iu±) 
(1, 1, 1, 1, 1; u, u, u, u, u) 
(u, u, u, u, u; — 1, — 1, — 1, — 1, — 1) 


Cj = cos(j7r/5) it = G+{a,b) 
Sj — sin(jV/5) v± = G±{—cia,C2b) 
F(a,(3) = y/5(a + /3) 2 - ±a(5 w± = G±(c 2 a,-ci6) 
E±(a,(3) = a + /3±F(a,/3) G±(a, (3) = [a - p ± F(a, /3)]/2(a + 0) 



TABLE II: Eigenvalues and eigenvectors of 5 qubits arranged on the vertices of a regular pentagon. 



n 


V 


Eigenvalues 


Eigenvectors 


1 


1,5 
2,4 
3 


a — b — c 
—a — b + c 
-2a + 26- c 




(1,-1, -2,, -1,,1,2);(1, 1,0,-1,-1,0) 
(-1,-1, 2,-1,-1, 2);(1,-1,0, 1,-1,0) 
(-1,1,-1,1,-1,1) 




6 


2a + 2b + c 




(1,1,1,1,1,1) 


2 


1,5 
2,4 

3 

6 


±Vb 2 + 3a 2 

V 

-2b 
2b 


(1, -1, -2, -1, 1, 2; 0, -w±, -w±, 0, u±, u±; 0, 0, 0) 
(3,3,0,-3,-3,0; 2u± , tt± , — u± , — 2u± , — u± , u± ; 0, 0, 0) 
(-1,-1,2,-1,-1,2; -2x, x, x, -2x, x, x; -y, 2y, -y) 
(1, -1, 0, 1, -1, 0; 0, -x, x, 0, -x, x- -y, 0, y) 
(-1,1, -1,1, -1,1; 0,0, 0,0, 0,0; 0,0,0) 
(0,0, 0,0, 0,0; -1,1, -1,1, -1,1; 0,0,0) 
(1, 1, 1, 1, 1, l;r, r,r,r,r,r;s,s,s) 






v A + 2bv 1 + (4ac - b 2 - 3a 2 - \c 2 )v + 26(6 2 - a 2 + 4ac) =0 x = (v 2 +bv- 2b 2 )/[v(a - 2c) - 2ab\ 
fj, 3 -46^ 2 -4(2ac + 6 2 + 3a 2 + c 2 )^ + 166(6 2 - a 2 - 2ac) =0 y = 2(a^ + 26c)/Ka - 2c) - 2a6] 
r = (^ 2 - 26^ - 86 2 )/2[^(a + c) + 4a6] s = 2(afi + 2bc) /[n{a + c) + 4ab] u± = -b±^Jb 2 + 3a 2 


3 


1,5 


a — 6 — c+(a + 6)m± 




(m±,-m±, -2m±, -m±,m±,2m±; 1, -1, -2, -1, 1,2; -1, -2, -1, 1,2, 1;0,0) 
(m±, m±, 0, — m±, — m±, 0; 1, 1, 0, —1, —1, 0; 1, 0, —1, —1, 0, 1; 0, 0) 




1,5 

2,4 


6 =p a =p c 




(0, 0, 0, 0, 0, 0; -2, T l, 1, ±2, 1, Tl; ±1, -1, T2, -1, ±1, 2; 0, 0) 




2,4 


c — a — o+(a — b)n± 




(0, 0, 0, 0, 0, 0; 0, -1, T l, 0, ±1, 1; 1, ±1, 0, T l, -1, 0; 0, 0) 
(n±, n±, -2n±, n±, n±, -2n±; 1, 1, -2, 1, 1, -2; 1, -2, 1, 1, -2, 1; 0, 0) 
(n±, -n±, 0, n±, -n±, 0; 1, -1, 0, 1, -1, 0; -1, 0, 1, -1, 0, 1; 0, 0) 




3 

3 
6 

6 


O- 

q=c±2a-26 




(-p-,p-,-p-,p-,-p-,p-;i,-i, 1,-1, i,-i; -1, 1,-1, 1,-1, i; -q-, q-) 
(o,o,o,o,o,0;tM,tM,tM;-i,±i,-i,±i,-i,±1;0,o) 

(p+,p+,p+,p+,p+,p+; 1, 1, 1, 1, 1, i; 1, 1, 1, 1, 1, i; <?+,?+) 






cr^ - (2b ± 3c ± 20)0-^ - [4(a 2 + b 2 ) + (2a ± 2b - c) 2 ]cr± ± 3c[c 2 + 2(±6c =F 4afe + ac)] = 
p± = 2[(26± a)a± + 3ca]/[a 2 : =F 2ccr ± - 3c 2 ] 
q± = 6[aa± =F ac + 2fec]/[a 2 =p 2ccr± - 3c 2 ] 
m± = [2c + 6-a± ^/(2c + 6-a) 2 + 8(a + b) 2 ]/[2(a + 6)] 
n± = [a + fe-2c± V(2c-6-a) 2 +8(a-b) 2 ]/[2(a-o)] 



TABLE III: Eigenvalues and eigenvectors of 6 qubits arranged on the vertices of a regular hexagon. 



N 


V 


Energy 


nearest neighbours 


linear dipoles, 


linear quadrupoles, 








only 


all neighbours 


all neighbours 


5 


1, 4 


.618a- 1.6186 


.618 


.239 


.474 




2, 3 


-1.618a + .6186 


-1.618 


-1.473 


-1.563 




5 


2a + 26 


2 


2.468 


2.178 


6 


1, 5 


a — 6 — c 


1 


.683 


.905 




2, 4 


—a — 6 + c 


-1 


-1.067 


-1.033 




3 


-2a + 26- c 


-2 


-1.741 


-1.902 




6 


2a + 26 + c 


2 


2.509 


2.159 



TABLE IV: n = 1 energies in the small sample limit (units of the static interaction energy between nearest neighbours). 
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iV 


AG (JV) 


n = 2 eigenvectors 


AG) \ 

(<") 


(in) 




(units of Viv) 




(units of Vat) 


(units of 7) 


2 


1 


(1) 





2 


3 


2 


(1.1,1) 


2 


4 


4 


2.354 


(1, 1, 1, l,x,x) 










x = 1.248 - .283i/Viv 


3.204 


5.930 






a; = -1.602 - .363i/V*r 


-2.496 


.070 


5 


2.472 












i- = 1.356 + M9i/V N 


3.823 


7.821 






x = -.737- .352i/Viv 


-1.351 


.179 


6 


2.511 


(1, 1, 1, 1, 1, l;x,x,x,x,x,x;y,y,y) 










x = 1.417 - 1.261i/V)v; y = 1.548 - 1.732i/Viv 


4.162 


9.687 






a; = -.114 - .373i/Viv; y = -1.087 - .550i/Vjv 


-.290 


.293 






x = -1.940 - .312i/V N ; y = 2.252 + .200i/V N 


-3.100 


.020 


7 


2.518 


(1, 1, 1, 1,1, 1, l;x,x,x,x,x,x,x;y,y,y,y,y,y,y) 










x = 1.45 - 1.602i/V)v; y = 1.66 - 2.738i/V N 


4.358 


11.534 






x = .314 - A33i/V N ; y = -.877 - .765i/V N 


.570 


.410 






x= -1.416- .379i/V N ; y= .634 + mii/V N 


-2.410 


.056 


8 


2.515 


(x,x,x,x,x,x,x,x; 1, 1, 1, 1, 1, 1, 1, l;y,y,y,y,y,y,y,y;z,z,z,z) 










a; = + 59.68i/VKr; y = 1.177 - .949i/V)v; z = 1.237 - lAUi/V N 


4.478 


13.37 






a; = 1.635 + 69.72i/Viv; y = - .0714i/V)v; z = -1.809 - 132.87i/Viv 


1.248 


.528 






x = -1.0914 - 285.8i/Vjv; y = -4523 + .0204i/Viv; z = -1.284 + .453i/Viv 


-1.660 


.097 






x = -.490 - 63.44i/Viv; y = -1.33 - 43.48i/Vjv; z = 1.451 - 34.9H/Vjv 


-3.322 


.008 


n 


z.oUo 


(1, 1, 1, 1, 1, 1, 1, 1, 1; x, x, x,x,x,x,x,x,x;y,y,y,y,y,y,y,y,y; z,z,z, z, z, z, z, z, z) 










x = 1.493 - 2.834i/Viv; y = 1-793 - 5.665i/V;v; z = 1.939 - 7Mli/V N 


4.559 


15.19 






x = .822 - .583i/V N ; y = -.144 - 1.184i/VW; z = -1.015 - 1.518i/V N 


1.780 


.644 






x = - .496 - A29i/V N ; y = -.916 - .396i/V N ; z = .713 - .0043i/Viv 


-.958 


.140 






x = -1.713 - .386i/VW; y = 1.538 + .177i/V N ; z = -.619 - .222i/Vjv 


-2.874 


.025 



TABLE V: Frequency shifts and widths of the exciton and biexciton energy levels (long- wavelength limit). 



N 


Frequency shifts 


Half- widths 


Relative Intensities 




(units of Vn) 


(units of 7) 




2 


-1 


4 


1 


3 





7 


1 


4 


0.85 


9.93 


.988 




-4.85 


4.07 


.012 


5 


1.351 


12.821 


.978 




-3.823 


5.179 


.022 


6 


1.651 


15.687 


.969 




-2.801 


6.293 


.029 




-5.611 


6.020 


.002 


7 


1.840 


18.534 


.961 




-1.948 


7.410 


.034 




-4.929 


7.056 


.005 


8 


1.964 


21.370 


.9550 




-1.266 


8.528 


.0377 




-4.174 


8.097 


.0069 




-5.836 


8.008 


.0006 


9 


2.052 


24.190 


.9494 




-.727 


9.644 


.0403 




-3.465 


9.140 


.0088 




-5.381 


9.025 


.0016 



TABLE VI: Biexciton excitation frequency shifts, natural line widths, and relative intensities (long- wavelength limit). 



